Combination of genetically diverse Pseudomonas phages enhances the cocktail efficiency against bacteria

Phage treatment has been used as an alternative to antibiotics since the early 1900s. However, bacteria may acquire phage resistance quickly, limiting the use of phage treatment. The combination of genetically diverse phages displaying distinct replication machinery in phage cocktails has therefore become a novel strategy to improve therapeutic outcomes. Here, we isolated and studied lytic phages (SPA01 and SPA05) that infect a wide range of clinical Pseudomonas aeruginosa isolates. These relatively small myophages have around 93 kbp genomes with no undesirable genes, have a 30-min latent period, and reproduce a relatively high number of progenies, ranging from 218 to 240 PFU per infected cell. Even though both phages lyse their hosts within 4 h, phage-resistant bacteria emerge during the treatment. Considering SPA01-resistant bacteria cross-resist phage SPA05 and vice versa, combining SPA01 and SPA05 for a cocktail would be ineffective. According to the decreased adsorption rate of the phages in the resistant isolates, one of the anti-phage mechanisms may occur through modification of phage receptors on the target cells. All resistant isolates, however, are susceptible to nucleus-forming jumbophages (PhiKZ and PhiPA3), which are genetically distinct from phages SPA01 and SPA05, suggesting that the jumbophages recognize a different receptor during phage entry. The combination of these phages with the jumbophage PhiKZ outperforms other tested combinations in terms of bactericidal activity and effectively suppresses the emergence of phage resistance. This finding reveals the effectiveness of the diverse phage-composed cocktail for reducing bacterial growth and prolonging the evolution of phage resistance.


Results and discussion
Phage morphological and biological properties of phages SPA01 and SPA05. The lytic phages named SPA01 and SPA05 were isolated from a canal in Bangkok, Thailand, using P. aeruginosa strain PAO1 as a parental bacterial host. Both phages were capable of forming 2-4 mm clear plaques on the PAO1 lawns (Fig. 1A,E), which the relatively large clear plaque appearance indicated the potent lytic activity of the phages against the bacteria. The morphology of phages SPA01 (Fig. 1B) and SPA05 (Fig. 1F), as examined by TEM and classified based on the criteria proposed by Ackermann 26 , revealed that both phages belong to the family Myoviridae due to the presence of an icosahedral head: 74 ± 9 nm and 75 ± 4 nm in diameter of SPA01 and SPA05 (n = 4), and a contractile tail: 163 ± 10 nm and 164 ± 5 nm in length of SPA01 and SPA05 (n = 4).
To investigate whether these phages are appropriate for medical and biotechnological purposes, their biological properties, including host range, adsorption rate, one-step growth curves, and phage tolerance were evaluated. The host range of phages against various Pseudomonas spp., including lab and clinical isolates, was first determined and the result showed that most of the clinical isolates were susceptible to phages SPA01 (10/13) and SPA05 (7/13) (Table 1). However, both phages were unable to infect the P. stutzeri, P. mendocina, P. fluorescens, and P. putida strains. This suggests that both phages are likely specific to a relatively wide spectrum among P. aeruginosa species, highlighting their potential use in medical treatment. To measure the infectivity of phages, the efficiency of plating (EOP) was performed. Of all strains tested, while none of them showed a higher EOP value than the indicator host strain, P. aeruginosa PAO1, the majority of EOP values obtained from phage SPA01 were higher than those obtained from phage SPA05 ( Table 1), indicating that SPA01 exhibited higher infectivity towards tested strains.
In aspects of adsorption, the majority of phages SPA01 (90%) and SPA05 (92%) adsorbed to the bacterial host cells within 10-15 min (Fig. 1C,G), followed by a 30-min latent interval for replication inside the cells (Fig. 1D,H). Both phages then burst the bacterial host, which resulted in approximately 218-240 virions per cell (Fig. 1D,H), highlighting the potency of these phages in therapeutic aspects since phages with high burst sizes are favorable [27][28][29] . The phage tolerance test showed that both phages were relatively stable over a wide range of temperatures and pHs ( Fig. S1A and S1B). The phages withstood the high temperature of 70 °C and did not lose much of their infectivity when exposed to acidic and basic conditions (pH 4-10), indicating the high viability of phages SPA01 and SPA05 as resilient biocontrol agents ( Fig. S1A and S1B). Altogether, these promising characteristics of phages SPA01 and SPA05 render them potential candidates for both therapeutics and biocontrol of P. aeruginosa.
Genomic characterization of phages SPA01 and SPA05. Apart from desirable phenotypic characteristics of phages that result in high infectivity against the targeted bacteria, genetic information is also crucial in selecting an appropriate, non-harmful phage for therapeutic use. To determine whether the undesirable genes associated with antibiotic resistance, bacterial virulence, and the lysogenic life cycle were present in the phage genomes, we conducted whole-genome sequencing, employed SPAdes 30 for de novo assembly, and annotated all open reading frames (ORFs) throughout the phage genomes. Genome analysis revealed that phages SPA01 and SPA05 contained 93,536 bp and 93,656 bp of double-stranded DNA genome, respectively (Fig. 1I,J); therefore, both phages were considered relatively small phages due to the genome size of less than 200 kb 31 . The %GC contents of the SPA01 and SPA05 genomes were 49.34% and 49.49%, respectively. The schematic genome maps of phages SPA01 and SPA05 created by the CG view server 32 showed the outer ring representing coding sequence locations (CDS) in blue and tRNAs in this ring in yellow. Genome annotation of phages using PHASTER 33 (Fig. S3). Moreover, the intergenomic similarities was below the threshold (95%) of identification of the similar species 50 . Therefore, these phages were classified as a novel phage in the genus "Pakpunavirus" (Fig. S3B). The ORFs were Figure 1. Morphological, biological, and genomic characteristics of phages SPA01 and SPA05. Plaque morphology of phages SPA01 (A) and SPA05 (E). Transmission electron micrographs of phages SPA01 (B) and SPA05 (F). Scale bar equals to 100 nm. Adsorption assays within 25 min of phages SPA01 (C) and SPA05 (G) with P. aeruginosa strain PAO1. One-step growth curve of phages SPA01 (D) and SPA05 (H) in P. aeruginosa strain PAO1 during a window of 70 min. Schematic whole genome maps of phages SPA01 (I) and SPA05 (J). The innermost circles colored in green and purple indicate the positive and negative GC skew, respectively. The open reading frames (ORFs) are indicated in blue color with arrows indicating the ORF direction. The functional annotation of these ORFs in phages SPA01 and SPA05 is shown in Table S1. The data shown in (C,D,G,H) represent the mean ± standard deviation of at least triplicates. www.nature.com/scientificreports/ mainly annotated as structural proteins, nucleotide metabolism and DNA replication-related proteins, and bacterial lysis proteins, as indicated in Table S1. Importantly, undesirable genes involved in bacterial virulence, toxins (exotoxin A, pyocyanin, hydrocyanic acid, enterotoxin etc.), or lysogeny were not identified within the phage genomes, assuring the lytic nature and implying safety for use of these two phages.
Potential of phages SPA01 and SPA05 in P. aeruginosa suppression and formation of phage-resistant bacteria. Since both phages SPA01 and SPA05 exhibited desirable characteristics of a good candidate for bacterial growth control, an in vitro lysis assay was performed to evaluate the lytic activity of phages SPA01 and SPA05 against P. aeruginosa at different multiplicity of infections (MOIs). The results showed that both phages displayed efficient killing activity against the bacteria, even at the very low MOI of 0.01 ( Fig. 2A,B). Although no bacterial regrowth was observed throughout a window of 240 min in both phages, viable bacteria could be detected at a longer timepoint ( Fig. 2A-D), suggesting the emergence of phage-resistant bacteria. SPA01 and SPA05 resistant P. aeruginosa isolates, in particular, were observed at 12 h and found to be approximately 2-log increased at 24 h in all tested MOIs (Fig. 2C,D). The revival frequency of bacteria toward phages SPA01 and SPA05 was quite comparable, which equals 0.0015-0.0017. This finding is unsurprising since it has been well documented that the use of phage, especially in vitro studies, normally leads to an early emergence of resistance [45][46][47] . Notably, despite resistance emergence during phage treatment, some studies have shed light on the possible evolutionary tradeoff of phage-resistant strains in P. aeruginosa that can lead to attenuation of bacterial virulence 46,48 or hypersensitivity to antibiotics 49,50 that would be beneficial for following treatments.
To further explore what possible resistance mechanisms the hosts exhibit against phages SPA01 and SPA05, we first isolated 10 P. aeruginosa strains that were resistant to each of the phages and examined the host adsorption rate, which is one of the most common anti-phage mechanisms 51 . Interestingly, the percentage of adsorption of each phage to all corresponding phage-resistant strains substantially dropped (Fig. 2E), suggesting that the mechanism of phage resistance in these strains possibly occurs via receptor modification. However, since phage resistance mechanisms can occur through the host immune system or genetic modification 9,52,53 , further investigation will be required to explore other possible resistance mechanisms in these isolates. We then investigated whether phage resistance to one phage affects susceptibility to another phage by measuring phage adsorption when cross-infected with the resistant isolates. The results demonstrated that all SPA01-resistant isolates were also resistant to phage SPA05 and vice versa, as indicated by the decreased phage adsorption percentage on the resistant strains (Fig. 2E). This finding suggested that these phages are closely related and infect the bacterial host through a similar adsorption mechanism. It has been shown that type IV pili (T4P) and lipopolysaccharide (LPS) O antigen are the common binding sites for P. aeruginosa phages [39][40][41][42]54,55 . Our VIRIDIC data revealed that phages SPA01 and SPA05 are classified as Pakpunavirus (Fig. S3). Therefore, they are assumed to recognize LPS Table 1. Bactericidal spectrum and efficiency of plating (EOP) of the phages. *Efficiency of plating (EOP) values were classified as highly productive (≥ 0.5), medium productive (0.1 ≤ EOP < 0.5), low productive (0.001 < EOP < 0.1) or inefficient (≤ 0.001). www.nature.com/scientificreports/ during the adsorption step, similar to other previous studied Pakpunavirus (Table S2). Thus, the impairment of phage adsorption in these resistant strains might affect the efficiency of the phage cocktail that is composed of SPA01 and SPA05. However, further investigation into the receptor, which serves a role in phage recognition, is needed to gain insights into the adsorption step of phages SPA01 and SPA05. It is well known that closely related phages are not compatible candidates for being the major components in a phage cocktail due to the possibility of cross resistance, which could jeopardize the cocktail's overall effectiveness. In order to diminish the emergence of phage-resistant bacteria, one of the strategies for cocktail design is to mix phages that are genetically divergent as they are assumed to display different replication mechanisms, such as targeting different host range and hijacking different bacterial molecular targets, particularly in this study as an example, recognizing different host cell receptors. This idea was previously supported by a study demonstrating that a combination of 5 genetically diverse phages against Mycobacterium tuberculosis greatly reduced persistence and phage-resistant bacteria 56 . However, insights into the replication machinery of these 5 genetically diverse Mycobacterium phages has still been unexplored.
Phages SPA01 and SPA05 are genetically and mechanistically diverse from jumbophages PhiKZ and PhiPA3. With the previous evidence that genetically diverse phages might serve roles in prolonging the emergence of phage resistance in the phage cocktail 56 , we reasoned that nucleus-forming Pseudomonas jumbophages that are believed to display a complicated infection machinery and could recognize a different phage receptor on the host cells would be beneficial as a candidate for the cocktail with phages SPA01 and SPA05. We first confirmed whether they are genetically different and investigated how they are phylogenetically clustered. The whole genome alignment of the phages demonstrated that the phage PhiKZ genome has a high degree of similarity in genome organization to that of phage PhiPA3, while the genomes of the jumbophages are organized differently from the genomes of the phages SPA01 and SPA05 (Fig. 3A). Phylogenetic analysis of the whole genome sequence further revealed that phages PhiKZ and PhiPA3 were closely related and grouped together in the clade of nucleus-forming jumbophages and speciated from other small phages, including phages SPA01 and SPA05 (Fig. 3B). Due to the large difference in genome size and organization among the phages, as www.nature.com/scientificreports/ expected, this result confirmed that our phages SPA01 and SPA05 are genetically diverse from the jumbophages PhiKZ and PhiPA3. The viral replicating or hijacking machinery, as we called it "Mechanism of pre-killing (MOK)", can be studied through the bacterial morphology that is triggered during phage infection. Through the morphological change, we have investigated how nucleus-forming jumbophages replicate inside bacteria 57-60 and predicted the bacterial fundamental metabolism that was interfered by phages 61 . Since we speculated that the small phages SPA01 and SPA05 might exhibit a distinct replication process during infection differing from the jumbophages PhiKZ and PhiPA3, as corresponded to their genetic divergence, we performed a single cell infection assay to observe if the morphological changes as triggered during the infection by these phages are different. The result revealed that all phages induced morphological changes in bacteria when compared to uninfected cells (Fig. 3C). In particular, SPA01-and SPA05-infected cells contained two compact nucleoids, while PhiKZ-and PhiPA3-infected cells were slightly longer and contained only one condensed nucleoid close to midcell (Fig. 3C). The morphology of PhiKZ-and PhiPA3-infected cells was well in agreement with our previous studies that revealed the presence of the phage nucleus that is centrally located near midcell during jumbophage infection 62 . We next quantitated the cytological profiles of each group at the single-cell level using morphological parameters reported previously 63,64 to test if the morphological changes caused by each phage are separated or clustered together. The result showed that SPA01-and SPA05-infected cytological profiles are clustered together and separated from PhiKZ-and PhiPA3-infected cells (Fig. 3C), suggesting that the two groups of phages trigger different morphological changes of cells during infection, possibly through different MOKs. Further studies will be needed to gain better understanding regarding MOKs of phage SPA01 and SPA05, in which it will provide insights into correlation between the MOK and the apparent morphological change. In summary, this evidence agrees with www.nature.com/scientificreports/ the evolutionary divergence of the phages, thus reassuring that these phages (SPA01 and SPA05) and the nucleusforming jumbophages are divergent.
The combination of the nucleus-forming jumbophages with either phages SPA01 or SPA05 efficiently suppresses the growth of bacteria and reduces the frequency of bacterial revival. During replication, PhiKZ-like viruses exhibit an orchestrated infection machinery by encoding a tubulin homolog named PhuZ to organize a nucleus-like compartment at midcell and deliver procapsids to dock on the phage nucleus surface for encapsidation [57][58][59][60]65,66 . The phage nucleus serves a role in partitioning proteins according to function and in protecting the phage nucleic acids against broad bacterial DNA-targeting immune systems such as clustered regularly interspaced short palindromic repeats (CRISPR) systems and restriction-modification (R-M) systems [65][66][67] . Due to their advantages over the small phages, as mentioned here, together with the idea of combining them with the phages SPA01 and SPA05 to formulate an efficient cocktail, we first tested whether the nucleus-forming jumbophages could lyse the phage SPA01 or SPA05-resistant isolates. Strikingly, both giant phages were able to kill all isolates that were resistant to each small phages SPA01 and SPA05, and efficiently suppressed the bacterial growth for at least 8 h (Fig. 4A,B). Results of the EOP also showed the high efficiency of the giant phages to lyse the resistant isolates in addition to their host bacteria (Fig. 4C). Seventeen phage-resistant isolates produced high levels of phage PhiPA3 production, having EOPs greater than 0.5, while thirteen strains could be classified as having high phage PhiKZ production (Fig. 4C). These results suggested that the giant phages (PhiKZ and PhiPA3) recognize different receptors on the target cells from those small phages (SPA01 and SPA05), making them a potential candidate for the phage cocktail formula. Due to the therapeutic potential of giant phages PhiKZ and PhiPA3 as they kill various clinical P. aeruginosa strains ( Table 1, Fig. S2) and are capable of infecting the phage SPA01 and SPA05-resistant isolates ( Fig. 4A-C), various combinations between the giant phages and the small phages (SPA01 and SPA05) were then formulated to make phage cocktails. The cocktail mixtures included 6 formulas as follows: (1) SPA01-SPA05, (2) SPA01-PhiKZ, (3) SPA01-PhiPA3, (4) SPA05-PhiKZ, (5) SPA05-PhiPA3, and (6) PhiKZ-PhiPA3. The cultures of P. aeruginosa PAO1 were grown in the presence of the phage cocktails at MOI 1 and the viable bacterial number was counted at 24 h and 48 h. The result showed that the remaining bacterial number that survived through the phage treatment was detectable in all formulas, suggesting the presence of phage resistance isolates in all conditions. However, the number of viable bacteria when treated with the small phage cocktail (SPA01-SPA05) was around 4 to 6 -log higher than those when treated with the combinations that contained a giant phage. The SPA01-SPA05 resistant bacteria outgrew faster than the others, with an around 2-log increase from 24 to 48 h (Fig. 4D). This result suggested that the closely related phages are not an appropriate combination in the cocktail since resistance to one might result in cross-resistance to another. Interesting, adding either one of the giant phages into the cocktail could efficiently suppress the regrowth of bacteria, as the viable bacterial numbers detected at both time points were significantly reduced when compared to the small phage cocktail. Furthermore, adding the giant phages, particularly PhiKZ, to the combination completely suppressed the bacterial regrowth, resulting in a comparable bacterial number at 48 h to 24 h, except for the cocktail SPA05-PhiPA3 (Fig. 4D).
The revival frequency of P. aeruginosa PAO1 treated with each phage cocktail was then investigated to evaluate the durability of phage therapy. The results showed that the frequency of revival developed in the small phage cocktail significantly increased from 24 to 48 h. The increase in revival frequency was also observed in the cocktails containing PhiPA3, whereas the increased revival was not observed in the cocktails that included PhiKZ (Fig. 4E), indicating that the giant phage PhiKZ is an appropriate ingredient for the cocktail with the small phages. Additionally, the cocktail of giant phages (PhiKZ-PhiPA3) alone also gave a promising outcome as it suppressed the bacterial regrowth from 24 to 48 h with unchanged revival frequency at 48 h when compared to the small phage cocktail (Fig. 4D,E). However, it is worth noting that, due to the survival of bacteria during the treatment of the PhiKZ-PhiPA3 combination that was significantly higher than that of the combinations of small and giant phages (Fig. 4D), this might be a result of co-infection of these closely related jumbophages in bacteria that might impair their infectivity as they interfered with each other through virogenesis incompatibility 68 .
It has long been known that phage cocktails consisting of different or distantly related phages are more effective in inhibiting the growth of targeted bacteria [69][70][71][72] . This is, in part, because different phages usually possess different mechanisms during replication 11 ; therefore, multiple resistances are required for the bacteria to fend off the infections caused by various phages in the cocktail. Thus, as observed here, combining the diverse phages that recognize different receptors on the target cells during phage entry in the cocktail can enhance the efficacy of the cocktail and decelerate the emergence of phage resistance. Implications for phage cocktail design. Conventional criteria for choosing phage cocktails are mainly based on empirical testing of different combinations of phages. Despite the fact that various methods for accelerating the process, such as high-throughput phage combination studies 73 , laborious work are still required for studying individual phage biology in order to properly select the phages for cocktail. Recently, the rise of genomic research and bioinformatic tools facilitates biologists in dissecting phage biology based on DNA sequence, which serve as a navigator toward choosing the right mix of diverse phages into cocktails 74 . However, the selection of phages for the cocktail is still largely based on genomic data. Here, we propose a complementary perspective for phage cocktail formulation via MOK-guided phage selection to choose genetically diverse phages for rapid cocktail design without the need for phage genomic data. Based on the morphological changes of bacterial cells triggered during phage infection together with the quantitation of the cytological profiles at a single-cell level, this information will provide additional criteria to the current phage selection checklist.
Although these MOK-based approaches through morphological changes would be beneficial in the selection of diverse phages for the rapid design of effective phage cocktails, there are some limitations in this study worth  www.nature.com/scientificreports/

Conclusion
We successfully isolated 2 virulent Pseudomonas myophages, SPA01 and SPA05, that exhibited potential therapeutic characteristics against various clinical strains of P. aeruginosa. Even though the phages efficiently suppressed the bacterial regrowth within 4 h, the phage SPA01 and SPA05-resistant bacteria subsequently emerged, possibly due to modification of a receptor where the phages bind during adsorption. The nucleus-forming jumbophages (PhiKZ and PhiPA3) that are distantly related to the phages SPA01 and SPA05 were thus selected as a potential candidate for formulating the phage cocktails. Through various cocktail formulas, we demonstrated that combining the jumbophages PhiKZ with either of the small phages improved the overall efficiency of the cocktail by better suppressing bacterial regrowth and delaying the emergence of phage resistance.

Materials and methods
Bacterial strains and growth conditions. Twenty-three Pseudomonas strains were employed in this study including 17 P. aeruginosa strains, 2 P. stutzeri strains, 2 P. putida strains, P. mendocina ATCC25411, P. fluorescens ATCC13525. A complete list of Pseudomonas strains can be found in Table 1. P. aeruginosa strains PSU-PA01-PSU-PA14, isolated from patients in internal medicine, Songklanagarind hospital, Thailand, were a generous gift from the RP lab. All bacteria were grown on LB broth (10 g/l tryptone, 10 g/l NaCl, and 5 g/l yeast extract) at 37°C followed by shaking at 250 rpm. When plated, cells were grown on LB containing 1.5% agar and incubated inverted at 37°C. This work has been reviewed and approved by Chulalongkorn University-Institutional Biosafety Committee ( Phage isolation, amplification, host range, and transmission electron microscopy. Water samples were collected from Khlong Samrong Bangkok, Thailand. Briefly, filtered water sample through a 0.45 µm membrane was mixed with 2X LB medium and then added log phase culture of P. aeruginosa PAO1 followed by incubation at 37 °C with shanking at 250 rpm for 24 h. After centrifugation, the supernatant was mixed with LB containing 0.5% agar powder and subsequently combined with the host. The mixture was overlaid onto LB plate and incubated at 37 • C for 24 h to observe plaque formation. The isolated plaque was collected and re-plated five times to ensure single phage purification. The purified phages were designated as phages SPA01 and SPA05. Phage stocks were maintained in SM buffer with gelatin (100 mM NaCl, 8 mM MgSO 4 ·7H 2 O, 50 mM Tris-Cl, pH 7.5 caution 0.01% (w/v) gelatin) and subsequently stored at 4 • C.
For phage propagation, 20 ml of the host culture (1 × 10 6 CFU/ml) was infected at a multiplicity of infection (MOI) of 1 and incubated for 24 h at 37 • C. The host cell was removed by centrifugation. The supernatant was filtrated (0.45 µm filter) and subsequently stored at 4 • C for further experiments. To determine phage titration, a tenfold serial of phage stocks was mixed with the host culture and poured into an LB soft agar plate. The plate was incubated overnight at 37 °C. The number of phages was enumerated by plaque observation and expressed as plaque forming unit (PFU) per ml.
Host range was identified by 23 Pseudomonas strains. The bacterial lawn was prepared with the bacteria and then 5 µl of phages were dropped onto the lawn of each host. The plate was incubated at 37 • C for 24 h to observe the lytic zone. The experiment was undertaken in triplicates. The efficiency of plating (EOP) values is calculated as the ratio of the PFU value of phage with susceptible bacterial strain and the phage with indicator bacterial strain.
Transmission electron microscopy (TEM) was undertaken to visualize phage morphology by negative staining method 75 . Briefly, fresh plaques were collected and suspended in SM buffer. The supernatant was placed on copper grids and negatively stained with 2% (w/v) uranyl acetate (pH 4.5). Phage morphology was observed with a Hitachi HT7700 transmission electron microscope at 80 kV.
Phage adsorption and one-step growth curve. The phage adsorption curve was determined by growing host cells to log phase (~ 1 × 10 8 CFU/ml), then infected with phages at an MOI of 0.01, and incubated at 37 °C. Aliquots of 1 ml were taken at an interval of every 5 min for 25 min and filtrated through 0.45 µm. The titer of unabsorbed phages was calculated from filtrates.
One-step phage growth curve was performed. Five milliliters of log phase host culture were centrifuged and resuspended in 1 ml of LB. Phages were transferred at an MOI of 0.01 and incubated at room temperature for 15 min for phage adsorption. The pellet containing infected cells was harvested and resuspended in 20 ml of LB broth and then incubated at 37 °C. Samples were taken at 10 min intervals for 70 min and subsequently filtrated through 0.45 µm. The titer of phages was immediately determined using the double-layer technique. This experiment was undertaken in triplicate. The latent period, burst time, and burst size of the phage were estimated as described elsewhere 76 . Effects of temperature, and pH on phage viability. Stability of phage at different conditions in SM buffer was determined. The various temperature (25 • C, 30 • C, 37 • C, 40 • C, 50 • C, 60 • C and 70 • C) and pHs (2, 4, 6, 8, 10, and 12) were prepared by diluting phage stock to a final concentration of 10 6 PFU/ml in 1 ml of SM buffer for 2 h. Subsequently, phage titers were measured as described above. The experiments were carried out in triplicate. To screen phage-resistant strains, P. aeruginosa PAO1 was grown at 37 • C in LB broth. When the bacterial culture reached OD 600 of 0.4, 10 µl of culture was mixed with 90 µl and 10 10 pfu/ml of each phage. The phageinfected culture was spread to LB agar plate. Then, plates were grown at 37 • C for 24 h. Ten randomly selected isolates were passaged five times on LB agar plate to confirm the stability of genetic changes. Resistance to SPA01 and SPA05 was confirmed in isolated clones by spot-test phage typing using phage suspension. Furthermore, the resistant strains were also tested the same way in terms of the susceptibility and the EOP to giant phages phiKZ and phiPA3. All tests were performed in triplicate.
Single-cell infection assay and image data analysis. The PAO1 was cultured to reach OD 600 ~ 0.4 and infected with phages at MOI of 5 followed by incubation for 20 min at 37 °C. The infected cells were fixed at a final concentration of 4% paraformaldehyde and incubated at room temperature for 15 min as described by Chaikeeratisak et al 57 . The fixed cells were harvested by centrifugation and then washed with 1 ml of 1 × PBS two times. The cells were resuspended in 1 × PBS and strained with fluorescent dyes (2 μg/ml FM 4-64 and 2 μg/ml DAPI). After centrifugation, 3 μl of cell suspension was loaded onto an agarose pad (1.2% agarose in 20% LB). The samples were visualized under Delta Vision Ultra High-Resolution Microscope.
Raw images from the fluorescent microscope were pre-processed on ImageJ software 77 . Then all individual cells were extracted morphological features by automatic cell analysis using Cellprofiler 4.0 software 78 . The morphological features set were selected based on the previous study 63 . Statistical calculation and Machine learning are carried out using scikit-learn library 79 in python. Briefly, cell profile data were transformed and normalized with Cube Root Transformation and Z-scored normalization 80 , respectively. Finally, the dimension of the dataset was reduced and visualized with Uniform Manifold Approximation and Projection (UMAP) 81 . To compare the uninfected and infected cell profiles, the uninfected cell of the phage treatment group was excluded by selecting the infected cell cluster on the UMAP plot that didn't overlap with the uninfected cell cluster.
Determination of bacterial revival frequency. The revival frequency of the bacteria was determined by a previous report with modification 82 . The PAO1 suspension was co-cultured with single phage or phage cocktail overnight at 37 °C with shaking at 150 rpm. The bacterial control without phages was incubated exactly as the infected cells. The number of bacteria was calculated based on the result of colony counting. The revival frequency was calculated using mean numbers of mutants in test samples divided by mean total numbers of control samples. Each independent experiment was repeated three times.
Phage genome sequencing and analysis. The phage DNA were extracted using a previous report with modification 83 . Briefly, purified phage was treated with DNaseI (2 mg/ml) and RNaseA (100 mg/ml) overnight at 37 • C. Subsequently, the lysis buffer (1 M Tris, pH 8.0, 0.5 M EDTA, pH 8.0, 10% SDS, and 10 mg/ml proteinase K) were added and incubated at 60 • C for 1 h. The phage DNA was extracted twice with phenol: chloroform: isoamyl alcohol (25:24:1) and chloroform. After centrifugation, the supernatant containing DNA was collected and precipitated by mixing of a 0.3 volume of 3 M NaOA and 1 volume of isopropanol followed by incubation for 2 h at − 20 • C. After the pellet was washed with 70% ethanol followed by centrifugation. The pellet was air-dried, dissolved with sterile distilled water, and then kept at −20 • C.The DNA concentration was estimated based on absorbance 260 by the nanodrop. The absorbance 260/280 ratio and agarose gel electrophoresis were performed to evaluate DNA quality and integrity.
The whole-genome sequencing of phage DNA was performed using short-read sequencing through MGISEQ-2000 (Beijing Genomics Institute, Beijing, China). The quality of reads was checked with FASTQC 84 and filtered using Trimmomatic 0.39 85 with default parameters. Spades 3.11.1 30 was used to assemble all reads into contigs. The open reading frames (ORFs) were identified using PHASTER 33 and GeneMarkS 34 following by annotation with Blastp of NCBI server. Antibiotic genes and virulence factors were identifed using ResFinder and VirulenceFinder 86 . The genome map of phages was drawn using the CG viewer server (http:// cgview. ca/) 32 . The genomic comparison of phages with closely related Pseudomonas phages was illustrated in the form of a linear figure using Easyfig application (http:// mjsull. github. io/ Easyfi g/ fles. html) 87 . Furthermore, tRNAs were predicted tRNAscan-SE (http:// lowel ab. ucsc. edu/ tRNAs can-SE/) 36 . The whole genome sequences were aligned in MAFFT online service 88 . The phylogenetic tree was generated in Geneious software (version 2022.2.2) using neighbour-joining method with 1000 bootstrap replications and visualized with Interactive Tree Of Life (iTOL) 89